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Early and sustained altered expression of 
aging-related genes in young 3xTg-AD mice 

V Gatta 1 ' 2 ' 7 , M D' Aurora 1 ' 3 ' 7 , A Granzotto 4 , L Stuppia 1 ' 2 and SL Sensi*' 3 ' 4 ' 5 ' 6 

Alzheimer's disease (AD) is a multifactorial neurological condition associated with a genetic profile that is still not completely 
understood. In this study, using a whole gene microarray approach, we investigated age-dependent gene expression profile 
changes occurring in the hippocampus of young and old transgenic AD (3xTg-AD) and wild-type (WT) mice. The aim of the study 
was to assess similarities between aging- and AD-related modifications of gene expression and investigate possible interactions 
between the two processes. Global gene expression profiles of hippocampal tissue obtained from 3xTg-AD and WT mice at 3 and 
12 months of age (m.o.a.) were analyzed by hierarchical clustering. Interaction among transcripts was then studied with the 
Ingenuity Pathway Analysis (IPA) software, a tool that discloses functional networks and/or pathways associated with sets of 
specific genes of interest. Cluster analysis revealed the selective presence of hundreds of upregulated and downregulated 
transcripts. Functional analysis showed transcript involvement mainly in neuronal death and autophagy, mitochondrial 
functioning, intracellular calcium homeostasis, inflammatory response, dendritic spine formation, modulation of synaptic 
functioning, and cognitive decline. Thus, overexpression of AD-related genes (such as mutant APP, PS1, and hyper- 
phosphorylated tau, the three genes that characterize our model) appears to favor modifications of additional genes that are 
involved in AD development and progression. The study also showed overlapping changes in 3xTg-AD at 3 m.o.a. and WT mice at 
12 m.o.a., thereby suggesting altered expression of aging-related genes that occurs earlier in 3xTg-AD mice. 
Cell Death and Disease (2014) 5, e1054; doi:10.1038/cddis.2014.11; published online 13 February 2014 
Subject Category: Neuroscience 



Alzheimer's disease (AD) is the most common form of 
dementia in the elderly population. 1 According to the 'amyloid 
cascade' hypothesis, amyloid dysmetabolism and formation 
of tau-containing neurofibrillary tangles (NFTs) are key steps 
in the disease. 2 Aging is also a critical factor for sporadic AD 
and a lively debate on whether the disease is actually driven 
by aging is animating the field. 3 The 'age-based hypothesis' 
postulates that initiating injury factors like head trauma, 
infections, vascular alterations, diabetes, or even genetic 
mutations, when operating on the aging brain, promote 
neuroinflammation and initiate the AD-related pathogenic 
cascade. 2 ' 3 Thus, brain aging along with neuroinflammation 
are crucial conditions on which pro-AD events must work to 
successfully initiate the disease. 

Genetic aspects of AD have been intensively investigated 
for decades (see Tanzi 4 for an extensive review on the topic); 
however, crucial information is still missing to successfully 
compose the puzzle. Microarray technology allows the 
simultaneous analysis of thousands of transcripts in a single 
experiment and is a useful approach for the investigation of a 
wide range of gene-related diseases. Gene expression 



studies in transgenic AD models have helped to unravel 
genetic factors influencing disease progression at well- 
defined stages as well as their relation with the development 
of cognitive deficits. 5 

Thus, in the present study we employed a whole genome 
microarray approach to investigate age-dependent gene 
expression profile changes in hippocampi obtained from 
young (3 months of age (m.o.a.)) and old (12 m.o.a.) 
3xTg-AD mice. This model offers the selective advantage of 
combining both amyloid (Aft)- and tau-dependent pathology 
and is a largely investigated and more comprehensive 
preclinical model of AD. 6 To control for age-dependent 
transcript modifications, we also investigated changes 
occurring in age-matched wild-type (WT) mice. 

In designing our experiment, gene expression of 12 m.o. 
3xTg-AD mice versus 3 m.o. 3xTg-AD mice were not 
considered as such design would have produced results 
indicative of changes that can be indistinguishably depending 
on a mix of two factors (aging and the AD-like background). 

The major aim of the study was to provide a better 
understanding of whether common genetic mechanisms are 
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Figure 1 Unsupervised hierarchical clustering analysis. Transcripts that are 
clustered according to their expression values (log ratios) are shown. Each row 
indicates a transcript. The nine columns depict three replicates for each of the three 
experimental conditions listed on the top of the figure. Quantitative changes in gene 
expression are shown in colors. Red and green indicate upregulated and 
downregulated transcripts, respectively. Black indicates no changes in expression. 
Missing data points are shown in gray, (a) Cluster of 35 transcripts that were found 
downregulated in 3 m.o. 3xTg-AD versus 3 m.o. WT mice and 12 m.o. 3xTg-AD 
versus 12 m.o. WT mice. The same genes were instead upregulated in 12 m.o. WT 
versus 3 m.o. WT mice, (b) Cluster of 59 transcripts that were found to be specifically 
upregulated in 3 m.o. 3xTg-AD versus 3 m.o. WT mice and in 12 m.o. WT versus 3 
m.o. WT mice. The same genes were instead downregulated in 12 m.o. 3xTg-AD 
versus 12 m.o. WT mice, (c) Cluster of 193 transcripts that were upregulated in all 
the experimental conditions (3 m.o. 3xTg-AD versus 3 m.o. WT mice and 12 m.o. 
3xTg-AD versus 12 m.o. WT mice as well as in 12 m.o. WT versus 3 m.o. WT mice), 
(d) Cluster of 76 transcripts that were found downregulated in 3 m.o. 3xTg-AD 
versus 3 m.o. WT mice and in 1 2 m.o. WT versus 3 m.o. WT mice. The same genes 
were instead upregulated in 12 m.o. 3xTg-AD versus 12 m.o. WT mice 



shared by aging and an AD-like background and whether 
overlapping pathogenic pathways can be identified in the two 
conditions. 

Results 

Cluster analysis of the investigated profiles revealed four 
different gene clusters (Supplementary Table 1 ), described in 
the following sections. 

Cluster A. Cluster A consists of 35 transcripts that, 
compared with age-matched WT mice, were found down- 
regulated in both 3 m.o. and 12 m.o. 3xTg-AD mice. The 
same gene set was found to be upregulated in 12 m.o. WT 
mice when compared with 3 m.o. WT (Figure 1a). The main 
biological functions associated with these genes are shown 
in Figure 2a and Table 1. 

Cluster B. Cluster B consists of 59 transcripts that, 
compared with 3 m.o. WT mice, were upregulated in 3 m.o. 
3xTg-AD mice and 12 m.o. WT mice. The same set was 
downregulated in 12 m.o. 3xTg-AD mice when compared 
with 12 m.o. WT mice (Figure 1b). These genes are 
implicated in key biological functions depicted in Figure 2b 
and Table 2. 

Cluster C. Cluster C consists of 193 transcripts that, 
compared with age-matched WT mice, were upregulated in 
3 m.o. and 12 m.o. 3xTg-AD mice. The same set was 
upregulated in 12 m.o. WT mice compared with 3 m.o. WT 
mice (Figure 1c). Functional analysis revealed that the 
overexpressed transcripts are associated with functions 
depicted in Figure 2c and Table 3. 

Cluster D. Cluster D consists of 76 transcripts that, 
compared with 3 m.o. WT mice, were downregulated in 
3 m.o. 3xTg-AD mice and 12 m.o. WT mice. These 
transcripts were upregulated in 12 m.o. 3xTg-AD mice when 
compared with 3 m.o. WT mice (Figure 1d). The main 
biological functions associated with these genes are provided 
in Figure 2d and Table 4. 

TaqMan qRT-PCR: microarray data validation. In order to 
validate the microarray results, quantitative real-time PCR 
(qRT-PCR) analysis was performed on RNA extracted from 
the same hippocampal samples employed for microarray 
experiments. Analysis was performed on three upregulated 
genes (BECN1, CST3, and GABRA5) belonging to cluster C. 
In all cases, qRT-PCR confirmed the microarray results 
(Figure 3). Moreover, expression of the employed house- 
keeping gene (GAPDH) was stable among all the samples. 

Interestingly, with the limitation and interpretative caution 
dictated by the low number of genes investigated by qRT-PCR 
(n = 3), we observed that young 3xTg-AD animals did already 
reach maximal expression that was not further changed, 
in a statistically significant manner, by aging (see 12 m.o. 
3xTg-AD animals). On the contrary, in WT mice, these 
changes were occurring in an age-dependent way. Interest- 
ingly, young 3xTg-AD mice actually showed expression 
values that were higher compared with old WT mice. 
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Figure 2 Biological functions as indicated by Ingenuity Pathway Analysis (IPA). Bar charts show results of IPA and indicate key biological functions modulated by genes 
selected in the four clusters described in Figure 1. (a) Cluster A, (b) cluster B, (c) cluster C, and (d) cluster D 



Discussion 

The main purpose of the study was to compare gene effects 
promoted by the pro-AD environment offered by the 3xTg-AD 
mice with those triggered, in WT mice, by senescence. 



The most relevant finding is the identification of profile 
modifications detected in young 3xTg-AD mice that closely 
resemble those occurring upon physiological aging in WT 
animals, thereby suggesting that 3 m.o. 3xTg-AD mice 
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Table 1 I PA functional analysis of cluster A genes 



Category 

Cellular movement 

Inflammatory response 

Cellular function and maintenance 

Cell cycle 

Lipid metabolism 

Molecular transport 

Cell death 

Tissue morphology 

Cell morphology 

Nervous system development and function 



P-value 

1 ,69E-03-4,63E-02 
1,69E-03-1,69E-03 
2,84E-03-3,66E-02 
3,38E-03-8,43E-03 
3,38E-03-3,38E-03 
3,38E-03-3,66E-02 
1,01E-02-3,84E-02 
1.18E-02-1.18E-02 
1,51E-02-4,92E-02 
3,66E-02-3,66E-02 



Molecules 

GPR182, PEBP1, CPLX3, TMSB 1 0/TMSB4X, HAS1 
TMSB 1 0/TMSB4X 

PEBP1, CPLX3, TMSB 1 0/TMSB4X, CDC42EP2, HAS1 

CKAP2 

PEBP1 

PEBP1, CPLX3 

TMSB 1 0/TMSB4X, CKAP2, DYSF, null 
HAS1 

PEBP1, TMSB 1 0/TMSB4X, ADCY6, VPS37C, HAS1 
CPLX3 



Table 2 I PA functional analysis of cluster B genes 



Category P-value Molecules 



Neurological disease 1 ,45E-04-4,55E-02 CTR9, NEUROD1, CD200, YWHAE, HSPA 1A/HSPA 1B, EGR1, F7, SSR1, SCOC, SNAP25, 

MFAP2, CDH2, SOX9, RSP02, PROCR, GABRA1 
Cell morphology 1 ,61 E-04-4,87E-02 NEUROD1, LDB3, YWHAE, PAM, HSPA 1A/HSPA 1B, EGR1, F7, OPHN1, SOX9, CDH2, 

NDUFS1, RSP02, PROCR, GABRA1 
Nervous system 1 ,14E-03-4,87E-02 NEUROD1, CD200, YWHAE, RIT2, HSPA 1A/HSPA IB, EGR1, ETV1, SNAP25, OPHN1, 

development and function CDH2, SOX9, GABRA 1, TACC1 

Cellular development 1 ,6E-03-4,43E-02 NEUROD1, SOX9, HSPA 1A/HSPA 1B, RIT2, EGR1, RSP02, F7, JAK3, SPRED2 
Tissue development 1 ,6E-03-4,75E-02 NEUROD1, MFAP2, SOX9,CDH2, CD200, EGR1, F7, JAK3, CERCAM 

Cell cycle 3,32E-03-3,64E-02 NEUROD1, SOX9, YWHAE, HSPA 1A/HSPA 1B, EGR1, JAK3 

Cell death 3,32E-03-4,95E-02 NEUROD1, CD200, YWHAE, HSPA 1A/HSPA 1B, EGR1, F7, DNAJB9, DSG1, NDUFS1, 

CDH2, SOX9, UBA7, PROCR, JAK3, TACC1 
Cell signaling 3,32E-03-3,59E-02 SOX9, CDH2, YWHAE, RIT2, IL10RB, F7, DOK3, SNAP25 

Lipid metabolism 3,32E-03-2,62E-02 Rdh1 (includes others,), F7 

Inflammatory response 1 ,65E-02-4,55E-02 CD200, EGR1, F7, DOK3, YTHDF2 



undergo premature modulation of aging-related genes 
(Figure 4). For clarity and to improve readability of the section, 
we are discussing results of the three experimental conditions 
as follows: 

• 3 m.o. 3xTg-AD versus 3 m.o. WT mice; 

• 12 m.o. 3xTg-AD versus 12 m.o. WT mice; 

• 12 m.o. WT versus 3 m.o. WT mice. 



3 m.o. 3xTg-AD versus 3 m.o. WT mice. Comparisons 
between expression profile changes in young 3xTg-AD mice 
versus 3 m.o. WT mice revealed overexpression of several 
genes that are depicted in clusters B and C (Figures 1 and 4). 
These genes modulate several mechanisms (discussed 
below) that are serving a pathogenic role in AD. 

[Ca 2 + ]j homeostasis. Deregulation of [Ca 2 + ]j is a key 
contributing factor for AD development and progression. 1 ' 7 
Perturbations in Ca 2+ handling by the endoplasmic reticulum 
(ER) and mitochondria as well as impaired intra/extracellular 
cation exchange contribute to the neuronal dysfunction and 
degeneration occurring in brain aging and AD. 7 In cluster C, 
we found five genes involved in [Ca 2 + ]j homeostasis 
(ATP2B1, CACNB4, RYR2, PKD2, and NDUFA9). ATP2B1 
encodes for the plasma membrane Ca 2 + -ATPase, a key 
system controlling high-capacitance/low-affinity Ca 2+ extru- 
sion, whereas CACNB4 encodes for the voltage-dependent 
L-type calcium channel subunit /?-4, a critical route for Ca 2 + 
influx. Altered expression and functionality of these proteins 



correlate with deregulation of [Ca 2 + ]j as well as with 
development of a seizure-prone phenotype. 8 CACNB4 is 
also associated with neuroprotective mechanisms in animal 
models of ischemic brain injury. 9 We also found changes in 
the expression of RYR2. RYR2 encodes for the ryanodine 
receptor 2, a major system controlling Ca 2+ release from the 
ER. RYR2 overexpression found in 3xTg-AD mice is in line 
with previous findings supporting an important role for 
the ryanodine receptor type 2 that is promoting altered 
ER-Ca 2+ release in presymptomatic 3xTg-AD mice. 10 This 
data set suggests a complex scenario in which enhanced 
[Ca 2 + ]j extrusion, mediated by ATP2B1, may compensate 
for impaired [Ca 2 + ]j buffering and excessive ER-Ca 2 + 
release operated by RYR2. 

Mitochondrial functioning. Alterations of mitochondrial 
functioning are crucial in physiological aging as well as in 
neurodegenerative conditions. 11 In clusters B and C, we 
found transcripts involved in the modulation of mitochondrial 
morphology and metabolism (TRAK1 and Ndufsl (cluster B); 
Hax1, NDUFA9, GDAP1, KIF1B, HADBH, PGAM1, Ppia, 
and DNM1L (cluster C)). GDAP1 and DNM1L encode 
for proteins that promote mitochondrial fission; PGAM1 
encodes for an enzyme of the glycolytic pathway that is 
involved in mitochondrial-mediated protective effects against 
A/?-induced cell death; TRAK1 plays a critical role in axonal 
transport of mitochondria. 12-15 Finally, Ndufsl encodes for 
the NADH-ubiquinone oxidoreductase 75 kD, a subunit that 
regulates the activity of mitochondrial complex I. The early 
and sustained overexpression of transcripts involved in 
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Table 3 I PA functional analysis of cluster C genes 



Category 



P-value 



Molecules 



Cell-to-cell signaling and 
interaction 

Cellular assembly and 
organization 



Cell morphology 



Organismal development 

Nervous system develop- 
ment and function 



Neurological disease 



4,94E-05-2,9E-02 
4,94E-05-2,9E-02 

6,66E-05-3,45E-02 

9,48E-05-2,9E-02 
1 ,75E-04-3,24E-02 

2,83E-04-3,34E-02 



Cell signaling 

Cellular compromise 
Cell death 



1 ,23E-03-3,34E-02 

1 ,23E-03-3E-02 
1 ,93E-03-3,45E-02 



Inflammatory response 5,57E-03-3E-02 



B2M, APBA2, VCAM1, GABRA5, USP14, JAK1, NRP2, MAPK1, KIF1B, CACNB4, CRK, HES1 
(includes EG: 15205), VDAC3, App, IRF1 (includes EG: 16362), SYBU, PMP22, DCC, RSU1 
B2M, APBA2, FHL1 (includes EG:14199), YWHAH, GPR12, MAPK1, KIF1B, HAX1, LPPR4, 
BOK, SNX3, CRK, BECN1, App, WASL, MYRIP, STARD13, DCC, MAP1LC3B, NRN1, 
VCAM1, NRP2, MED1 (includes EG:19014), RYR2, SEPT7, SMARCA5, YWHAZ, STK38L, 
MAPK9, GDAP1, VDAC3, TSG101, SYBU, PMP22, CETN3, Hmg1l1, HPRT1, DNM1L, ENC1 
B2M, APBA2, VCAM1, MAPK1, GPR12, NRP2, SEPT7, RYR2, STK38L, YWHAZ, CRK, 
GDAP1, HES1 (includes EG:15205), BECN1, VDAC3, App, TSG101, PMP22, Hmg1l1, 
STARD13, DCC, MAP1LC3B, DNM1L, NRN1 

PAFAH1B2, VCAM1, JAK1, MAPK1, NRP2, MED1 (includes EG:19014), PKD2 (includes 
EG:18764), MAPK9, BECN1, App, H2AFZ, TSG101, WASL, Hmg1l1 ,CST3,TIMP2 
B2M, APBA2, USP14, GABRA5, GPR12, MAPK1, YWHAH, KIF1B, ITM2B, LPPR4, CRK, 
HES1 (includes EG:15205), App, SOX2, CST3, DCC, BPNT1, NRN1, HEY1, VCAM1, NRP2, 
MED1 (includes EG:19014), CACNB4, STK38L, YWHAZ, MAPK9, VDAC3, SYBU, PMP22, 
null, Hmg1l1, HPRT1, DNM1L, ENC1 

GABRA5, YWHAH, MAPK1, ATP6V1C1, ATP2B1, KIF1B, LPPR4, SLMAP, SOX2, PDHA1, 
TTBK2, PGAM1, GABRB1, PKIA, PGRMC1, G3BP2, IMPA1, HIAT1, YWHAZ, SPARCL1, 
HPRT1, CYC1, TMEM66, RSU1, ORC4, ENC1, B2M, APBA2, CA2, FHL1 (includes 
EG:14199), ITM2B, PPP1CB, BECN1, App, UNC79, TBCEL, FARSB, HADHB, TMEFF1, 
PPP3CB, CST3, STARD13, DCC, null, JKAMP, VCAM1, NRP2, SLC2A13, CACNB4, RYR2, 
MAPK9, GDAP1, KIFAP3 (includes EG: 16579), HSPA12A, PMP22, AK5, null, DOCK8, PPIA, 
GLT25D2, DNM1L 

G3BP2, JAK1, MAPK1, ATP2B1, PKD2 (includes EG:18764), RYR2, CACNB4, STK38L, 
MAPK9, CRK, KIFAP3 (includes EG:16579), App, IRF1 (includes EG:16362), SOX2, CLK1, 
TRIM13, TTBK2, PPM1A, RSU1, TMEM9B, TIMP2 

VCAM1, MAPK1, SEPT7, PGAM1, CRK, GDAP1, DNM1L, TSG101, App, IRF1 (includes 
EG: 16362) 

PAFAH1B2, CA2, VCAM1, JAK1, MAPK1, KIF1B, MED1 (includes EG:19014), CUL1, YWHAZ, 
SMARCA5, MAPK9, BOK, HES1 (includes EG:15205), BECN1, KIFAP3 (includes EG:16579), 
TSG101, App, IRF1 (includes EG:16362), GNPTG, NT5C3, PMP22, AK5, Hmg1l1, DCC, 
PPM1A, TIMP2 

VCAM1, JAK1, Hmg1l1, PPIA, PGAM1, App, IRF1 (includes EG:16362), TIMP2 



Table 4 I PA functional analysis of cluster D genes 



Category 



P-value 



Molecules 



Cell-to-cell signaling and 
interaction 
Nervous system 
development and function 
Organ development 
Cellular development 

Cell cycle 

Inflammatory disease 
Lipid metabolism 
Cell death 

Cell morphology 
Cellular growth and 
proliferation 
Neurological disease 



9,74E-05-4,89E-02 

9,74E-05-4,99E-02 

3,38E-04-4,59E-02 
4,87E-04-4,78E-02 

8,76E-04-4,78E-02 

2,39E-03-3,47E-02 
3,14E-03-3,6E-02 
4,07E-03-3,21E-02 

4,07E-03-4,39E-02 
4,07E-03-4,99E-02 

8,12E-03-3,21E-02 



ST5, EPHB1, ARID 1 A, CCND3, BAIAP2, RGS4, STAT3 

ISLR2, NEUROD2, EPHB1, PRPF19, CCND3, ACTB, SWOB, BAIAP2, PDE4A, RGS4, 
STAT3, DI02, KLK8, RNF6, HIP1R 
EPHB1, CCND3, STAT3, DI02 

PRPF19, PA2G4, ELAVL3, PDE4A, RGS4, STAT3, DI02, BIN1, GPC1, NEUROD2, LSMD1, 
CCND3, SWOB, BAIAP2, BRD4, UBTF, TRRAP, SCAND1 

GPC1, MAP4, ARID1A, null, CCND3, PA2G4, EHMT2 (includes EG:10919), LAS1L, BRD4, 

STAT3, TRRAP 

CCND3, PDE4A, STAT3 

GPC1, MAP4, null, COTL1, HIP1R 

PRPF19, ARID1A, PA2G4, ACTB, PDE4A, RGS4, STAT3, CAPNW, KLK8, BIN1, ITPK1, 
GPC1, MAP4, LSMD1, CCND3, SWOB, VCP, UBTF, TEX261 

ISLR2, MAP4, ACTB, BAIAP2, VCP, SWOB, PDE4A, BRD4, RGS4, STAT3, TRRAP, RNF6 
PRPF19, PA2G4, ACTB, RGS4, STAT3, NAAW, BIN1, GPC1, LSMD1, CCND3, EIF4A1, 
SWOB, BRD4, PTPRS, LAS1L, UBTF, TRRAP 
SWOB, STAT3 



mitochondrial fission (GDAP1, and DNML1) and neuro- 
protection (PGAM1) is in line with recent findings indicating 
that mitochondrial fission is neuroprotective through reduc- 
tion of oxidative damages. 15 

Inflammatory response. Neuroinflammation is detrimental 
in brain aging and age-related cognitive decline. 16 The AD 
brain shows chronic, discrete, and microlocalized areas 
of inflammation as well as the appearance of activated 



microglia and reactive astrocytes in the vicinity of senile 
plaques. 1 Moreover, inflammatory-related genes and pro- 
teins (i.e., CD33, NLRP3) play a substantial role in A^-driven 
neurodegeneration and AD progression. 17,18 Cluster C 
revealed the presence of some transcripts (APP, IRF1, 
PPIA, VCAM1, TIMP2, HMGB1L1, PGAM1, and JAK1) that 
encode for proteins involved in the inflammatory response. 
It is noteworthy that in analogy with what was found for 
mitochondrial functioning, 3xTg-AD mice showed an early 
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activation of upstream (IRF1) and downstream (JAK1) 
proinflammatory transcripts. 

Additional genes involved in the inflammatory response 
(Hsp70, IL 10 receptor, and CD200) were also found (cluster D). 
However, these transcripts were downregulated (instead of 
upregulated when comparing 3 m.o. 3xTg-AD mice with 
3 m.o. WT mice). Among these, Hsp 70 and CD200are known 
to play an anti-inflammatory role. Thus, the overexpression of 
inflammatory genes appears to be coupled with simultaneous 
decrease in the expression of anti-inflammatory transcripts, 
thereby suggesting an overall proinflammatory status. 
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Figure 3 Validation of microarray gene expression data by qRT-PCR. Bar 
graph shows mRNA levels of CST3, GABRA5, and BECN1 measured by real-time 
PCR in young and old WT mice and 3xTg-AD mice as well as in old WT mice and 
3xTg-AD mice. Data are expressed as mean values of relative fold changes ± S.D. 
of three independent experiments performed in triplicate. The symbol '*' indicates 
that CST3, GABRA5, and BECN1 are significantly upregulated when compared 
with gene expression in 3 m.o. WT mice (P<0.05). The symbol 1# ' indicates that 
CST3, GABRA5, and BECN1 are significantly upregulated when compared with 
expression of these genes in 12 m.o. WT mice (P<0.01) 



Dendrite formation, synaptic plasticity, and memory. 

Age- and AD-related cognitive decline is associated with 
synaptic dysfunction. 19 In that respect, we found over- 
expression of several transcripts (CD200, Neurodl, SOX9, 
SNAP25, EGR1, OPHN1, ETV1, CDH2, GABRA1, and 
Spondl (cluster B); APBA2, APP, GABRA5, MAPK1, 
PMP22, SYBU, USP14, VDAC3, CST3, DCC, HES1, 
HEY1, HPRT1, SOX2, YWHAH, CRK, GPR12, HMGB1L1, 
STK38L, and VCAM1 (cluster C)) that are implicated in 
dendritic spine formation, modulation of synaptic functioning, 
and cognitive decline. Overexpression of transcripts encod- 
ing for proteins involved in synaptic activity and neuro- 
transmission may be interpreted as a protective event 
against AD-related synaptic dysfunctions. 

Overexpression of transcripts present in cluster C 
suggests a bimodal set of activities. One set of actions may 
counteract neuronal damage and promote protective effects 
aimed at maintaining cognitive abilities. For instance, APBA2 
is known to bind to APP and reduce A/? production, thereby 
counteracting memory deficits and promoting long-term 
potentiation (LTP) in Tg AD mice. 20 Similarly, CST3 up- 
regulation is protective against the neurotoxic effects of Aft 
oligomers and oxidative stress damage. This phenomenon 
has been reported in Tg AD models, in the brains of AD 
patients, and in those of elderly individuals. 21 The GABRA5 
transcript increase found in our 3xTg-AD mice lends 
support to the compensatory hypothesis as GABRA5 
has been reported to be downregulated with aging and this 
process has been linked to defective spatial memory 
performance. 22 

These protective actions are counteracted by the concur- 
rent upregulation of transcripts that, on the contrary, are 
associated with neuronal dysfunction. This is the case of APP 
and ERK2 transcripts. Given the aminoacidic composition of 
murine A/?, murine APP is marginally involved in the formation 
of toxic Aft species, but some studies have also proposed that 
the protein can favor ERK1,2 activation, tau phosphorylation, 
and the formation of NFTs. 23 Chronic ERK2 activation by A/? 
oligomers can be detrimental as early and sustained ERK2 
stimulation has been shown to promote cognitive decline in 
AD mice 24 and induce neuronal toxicity through caspase-3 
activation 25 
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Figure 4 Graphical synopsis of cluster analysis results. Synoptic view of all the transcriptome profiles that were found altered in the three experimental conditions. 
The main findings associated with observed gene expression changes are summarized. Note that, as depicted by the dark gray boxes, 3 m.o. 3xTg-AD versus 3 m.o. WT mice 
as well as 12 m.o. WT versus 3 m.o. WT mice show a similar pattern of expression changes in clusters B, C, and D. Also note that, only in the case of genes of cluster C, 
a common pattern of expression profiles is observed across all the three experimental conditions 
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It is interesting to note that in clusters A and D, we found 
downregulation of transcripts (PEBP1 and Complexin III 
(cluster A) and BIN, ACTB, BAIAP2, EphB, KLK8, ARD1, 
PDE4, Ehmt2, Neurod2, STAT3, SWOB, and CCND3 (cluster 
D)) that are involved in dendritic spine formation, modulation 
of synaptic functioning, and cognitive decline. In this respect, 
PEBP1 and Complexin III downregulation is associated 
with impaired cholinergic neurotransmission (PEBP1) 26 and 
synaptic transmission deficits (Complexin III), 27 the two 
mechanisms largely impaired in AD and in mild cognitive 
impairment. 28 

BAIAP2, EphB1, BIN1, and KLK8 are involved in synaptic 
plasticity. BAIAP2 (also known as IRSp53) modulates NMDA 
receptor-related synaptic transmission, LTP, as well as 
learning and memory. BAIAP2 knockout mice show impaired 
spatial learning and novel object recognition. 29 Thus, BAIAP2 
downregulation might be involved in the development of 
cognitive deficits affecting this model that has been also 
shown to suffer from significant decrease in BDNF levels. 30 
NMDAR deregulation has been associated with EphB 
(a family of ligand proteins) receptor deficits and linked to 
the appearance of neurological disorders including AD. 31 Still, 
in the set of synaptic plasticity-related genes that we found 
underexpressed, KLK8 seems important as the gene is 
involved in brain development, neuronal plasticity, and stress 
response. 32 Thus, KLK8 downregulation can favor aberrant 
CNS development. SWOB is a calcium-binding protein 
released by astroglial cells. Compelling evidence indicates 
that transgenic animals, overexpressing S100B, show 
neuronal loss and increased expression of the pro-apoptotic 
protein clusterin, a protein that is increased in hippocampi and 
frontal cortices of AD patients. 33 

In summary, our findings suggest that 3xTg-AD mice 
not only show early modulation of transcripts that exert 
detrimental effects on synaptic functioning but also develop 
changes that counteract this impairment. The process 
appears to go on chronically as the same pattern is present 
in 12 m.o. 3xTg-AD mice. 

Neuronal death, cell cycle, and autophagy. Aberrant 
expression of cell cycle- and autophagy-related molecules 
plays a pivotal role in the neuronal loss associated with AD 
and brain aging. 34 ' 35 Our analysis revealed overexpression 
of a group of transcripts (SPIN1, JAK3, Neurodl, SOX9, 
YWHAE, EGR1, and HSPA 1A/HSPA 1B (cluster B); Beclin-1, 
MAP1LC3-B, HADHB, AK5, CyC1, SIRT7 and ENC1 
(cluster C)) involved in autophagy, cell cycle, and cell death, 
thereby suggesting the potential activation of a complex 
network of neuronal death-related responses occurring in 
young 3xTg-AD mice. 

Early EGR1 upregulation has been reported in AD 
patients 36 and can contribute, along with other cell cycle- 
related genes, to NFT formation, a phenomenon that 
is also mediated by the increased expression of cell cycle 
proteins. 34 ' 37 Furthermore, overexpression of SPIN1 (a 
gene involved in cell cycle impairment) promotes chromo- 
some instability, a process leading to cellular senescence and 
apoptosis. 38 Sirt7, a less investigated member of the sirtuin 
family, is a nuclear protein expressed in the brain with a 
cellular distribution and functions that are still unknown. 39 



Experiments employing Sirt7 knockout mice or Sirt7- 
overexpressing cells have indicated the antiproliferative role 
of the protein and suggested that Sirt7 activity can improve 
tissue integrity in aging animals 40 

Beclin-1 and MAPI LC3-B are two autophagy-related genes 
whose overexpression might counteract pathology develop- 
ment in 3 m.o. 3xTg-AD mice. This hypothesis is consistent 
with recent findings indicating that early induction of 
autophagy reduces cognitive decline in 3xTg-AD mice 41 42 
Autophagy, in fact, clears out damaged proteins that are 
prone to aggregation, 43 and the process can be neuroprotec- 
tive by mediating the degradation of misfolded proteins. 44 
Intriguingly, in AD, autophagy can play a dual role. On one 
hand, it is protective by promoting the removal of A/? and tau 
aggregates and, on the other hand, autophagosomes are 
nevertheless sites of choice for accumulation and production 
of A£ 45 " 47 

Analysis of networks of genes involved in neuronal death 
also showed several transcripts that, compared with young 
WT mice, were found downregulated in young 3xTg-AD mice 
(ARID1A, BRD4, CCND3, EHMT2, GPC1, LAS1L, MAP4, 
PA2G4, PTGES3, STAT3, and TRRAP (cluster D)). Interest- 
ingly, aberrant expression of CCND3 has been reported to be 
closely associated with AD-related neuronal death occurring 
in specific hippocampal regions 48 

12 m.o. 3xTg-AD versus 12 m.o. WT mice. In order to 
evaluate how AD-like progression affects gene expression, 
we investigated profiles of 12 m.o. 3xTg-AD mice compared 
with 12 m.o. WT mice and matched this set of data with 
results of 3 m.o. 3xTg-AD mice that were compared with age- 
matched WT mice. The analysis revealed a similar over- 
expression of all transcripts of cluster C and downexpression 
of those of cluster A. Interesting differences appeared 
when analyzing transcripts belonging to clusters B and D 
(Figures 1 and 4). 

Cluster B transcripts, overexpressed in young 3xTg-AD 
mice, were downregulated in 12 m.o. 3xTg-AD mice when 
compared with age-matched WT mice, a finding that supports 
the hypothesis that aged 3xTg-AD mice may show deficits in 
counteracting disease progression. Underexpression of many 
transcripts found in cluster B has been previously reported in 
several models of aging and AD as well as in AD patients. This 
is the case for SNAP25 (encoding for a key protein involved in 
synaptic vesicle cycle and neurotransmission modulation), 49 
GABRA1 (encoding for a y-aminobutyric acid (GABA) 
receptor subunit downregulated in AD brains), 50 CDH2 
(encoding for N-cadherin, whose deregulation triggers 
apoptotic pathway activation), 51 and OPHN1 (encoding for a 
Rho-GTPase involved in controlling excitatory synapses 
maturation and plasticity). 52 

Cluster D transcripts, downexpressed in 3 m.o. 3xTg-AD, 
were instead overexpressed in aged 3xTg-AD. The upregula- 
tion of transcripts that are involved in the modulation of 
synaptic plasticity, neuronal proliferation, survival, differentia- 
tion, and regulation of cell division process might reflect a 
compensatory neuroprotective response (BAIAP2), although 
the net result suggests a more complex scenario (see, for 
instance, the role of KLK8). BAIAP2 downregulation is 
related to altered synaptic plasticity as observed in 
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BAIAP2-deficient mice, whereas its overexpression 
increases dendritic spine density and postsynaptic potentia- 
tion, ultimately enhancing cognition. 53 KLK8 deficiency 
promotes decreased synaptic plasticity, whereas, on the 
contrary, gain in KLK8 mRNA levels has been observed in AD 
brains and equally associated with detrimental effects on 
hippocampal functioning. 54 Inhibition of SWOB has positive 
effects on the A/? load and gliosis observed in AD patients, 55 
whereas SWOB overexpression has been associated with 
exacerbation of AD clinical signs. 56,57 Finally, altered expres- 
sion of cell cycle proteins (i.e., CCND3) in postmitotic neurons 
has been observed both in aging and in AD. 

Collectively, over- and underexpression of the afore- 
mentioned transcripts confirm a pattern of responses in 
which some genes exert detrimental effects that are in line 
with the development of AD-like pathology. At the same 
time, other transcripts appear to counteract these effects. 
Interestingly, this process appears to progress chronically as 
a similar pattern has been reported in 3 m.o. and in 12 m.o. 
3xTg-AD mice. 

12 m.o. WT versus 3 m.o. WT mice. Finally, we compared 
expression profiles of 12 m.o. WT mice with 3 m.o. WT mice 
to analyze the effects of physiological aging. 

Interestingly, changes driven by aging were largely over- 
lapping with those observed in 3 m.o. 3xTg-AD mice when 
compared with 3 m.o. WT mice. We, in fact, observed 
overexpression of several transcripts belonging to clusters B 
and C and downregulation of a gene data set present in 
cluster D (Figures 1 and 4). The functional correlates of these 
changes have been discussed above. 

Overall, these data support the hypothesis that the pro-AD 
environment offered by the 3xTg-AD mice accelerates and 
promotes early modifications of genes that are otherwise 
changed by aging in old WT mice. Moreover, comparison 
between 12 m.o. WT mice and age-matched 3xTg-AD mice 
showed expression changes that were similar only for genes 
belonging to cluster C. 

Interestingly, Ingenuity Pathway Analysis (IPA) analysis 
of this cluster revealed transcripts (HSPA12A, HIAT1, 
KIAA1409, NRP2, STARD13, UBB, and PKIA) that genome- 
wide association studies have recently shown to be involved 
in AD-related single-nucleotide polimorphisms (SNPs) as 
well as transcripts (APP, ATP6V1C1, BECN1, CST3, DCC, 
DNM1L, GABRA5, GABRB1, HIAT1, HSPA12A, KIAA1409, 
LPPR4, MAPK9, NRP2, PKIA, RYR2, STARD13, UBB, and 
YWHAZ) that are closely related to AD development. 

Finally, cluster A is the only one showing a unique behavior 
in 3xTg-AD mice, indicating changes that are specifically 
triggered by the pro-AD genetic background of the mice. It is 
intriguing to hypothesize that these genes have a functional 
role in AD progression and represent potential AD biomarkers. 
Among these, we observed a downregulation of Homer3 
and Dysferlin. Homer3 interacts with APP and inhibits A/? 
production. Thus, Homer3 downregulation can lead to 
enhanced A/? production. 58 Accumulation in senile plaques 
of Dysferlin, a protein involved in membrane repair, leads to 
defective neuronal repair. 59 Thus, it is conceivable that the 
Dysferlin downregulation that we found in the 3xTg-AD mice 
can exert a similar detrimental effect. 



Conclusions 

With the limitation of any Tg-AD model, our data support the 
hypothesis that aging and AD share some mechanistic 
similarities. In our Tg-AD mouse, the pro-AD environment 
promotes and accelerates changes that, upon aging, are 
otherwise occurring in WT animals. Although the 3xTg-AD 
mice harbor three mutations that do not occur, at the same 
time, in AD human brains, this cumulative pro-AD background 
modulates in Tg mice the early on expression of genes, the 
so-called BioAge genes (associated with neuronal loss, glial 
activation, and lipid dysmetabolism) that are in humans 
biomarkers for brain aging and found prematurely expressed 
in AD patients. 60 The same accelerated expression of aging- 
related genes in AD has been recently reported in a study that 
analyzed AD brains, 61 ' 62 thereby providing support for the 
heuristic value of 3xTg-AD mice as preclinical model of AD. 

Our data support and integrate the age-based hypothesis 
for AD, a concept that postulates that aging is required to set in 
motion the AD-related pathogenic cascade. Our findings 
suggest the presence of a feed-forward mechanism by which 
pro-AD factors can promote aging-related gene profile 
changes, thereby boosting and accelerating the disease 
process. Furthermore, this is the first study that offers an 
extensive analysis of the transcriptomic profile of the 3xTg-AD 
model, one of the most comprehensive preclinical models of 
the disease. The novel set of AD-related genes found in our 
analysis will, likely, further studies investigating the molecular 
roadmap of the disease and its modulation by aging. 

Materials and Methods 

Animal models. All the procedures involving animals were approved by the 
institutional ethics committee (CeSI protocol no.: AD-301) and performed 
according to institutional guidelines and in compliance with national (DL no. 
116, GU, Suppl. 40, 18 February 1992) and international laws and policies. The 
3xTg-AD (harboring AD- related human mutation APP (Swe), PS1 (M146V), and 
tau (P301L) transgenes) female mice (n = 2 at 3 m.o.a.; n = 2 at 12 m.o.a.) and 
age-matched female control WT mice (129SV x C57BL/6, n = 2 at 3 m.o.a.; 
n = 2 at 12 m.o.a.) were kept in a temperature-controlled room at 25°C under a 
12-h light/dark cycle and fed ad libitum with tap water and commercial chow. Mice 
at 3 or 12 m.o.a. were anesthetized and killed by decapitation, hippocampi 
excised, tissue transferred into RNA-later solution, and stored at - 80°C for RNA 
processing. 

Microarray analysis. Hippocampi were homogenized using a hand glass 
potter and total RNA extracted using the SVtotal RNA Isolation System kit 
following the manufacturer's instructions (Promega, Madison, Wl, USA). RNA 
purity and quantity were estimated using an Agilent 8453 spectrophotometer 
(Agilent, Santa Clara, CA, USA). RNA samples were not pooled. Using the Amino 
Allyl MessageAmp II aRNA Amplification Kit (Ambion, Grand Island, NY, USA) 
RNA was amplified and fluorescently labeled with Cy3-Cy5 cyanins following the 
manufacturer's instructions. Labeled RNA was hybridized on high-density array 
containing 31 802 mouse transcripts (Mouse OneArray Whole Genome DNA 
microarray v1, Biosense, Milan, Italy). In order to increase the experimental 
homogeneity, we performed three replicates with a dye swap and each experiment 
contained RNA from hippocampi of young (3 m.o.a.) and old (12 m.o.a.) 3xTg-AD 
mice versus age-matched WT hippocampal RNA as control for a total of nine 
experiments. After hybridization, Cy3-Cy5 fluorescent signals were detected with 
a confocal laser scanner 'ScanArray Express' (Packard Bioscience, Waltham, MA, 
USA) and analyzed with the 'ScanArray Express-MicroArray Analysis System 3.0' 
software (Perkin Elmer, Waltham, MA, USA). Raw data are stored in the GEO 
public database (accession number: GSE35210). Raw data were normalized using 
locally weighted scatterplot smoothing (LOWESS) algorithm. A transcript was 
considered differentially expressed when showing an absolute log-ratio value of 
> ±0.5, that is, 1.4 fold-change in transcript quantity. The transcript data set 
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was analyzed with gene clustering analysis using Cluster 3.0 (Stanford University 
Labs, Stanford, CA, USA). We evaluated only spots showing significant 
expression in at least 80% of the experiments. Data were centered to the mean 
of each mRNA and unsupervised average linkage hierarchical clustering of both 
the mRNAs and samples was performed using uncentered correlation metrics. 
Dendrograms and selected transcripts were visualized by employing TreeView 
software (Stanford University Labs). Gene biological functions were inferred using 
IPA software (Ingenuity Systems, Redwood City, CA, USA). 

Real-time PCR. Microarray results were validated by performing qRT-PCR 
analysis on RNA extracted from hippocampi obtained from the same samples 
used for microarray analysis. Gene expression of three genes found upregulated 
by microarray (BECN1, CST3, and GABRA5) was evaluated. The GAPDH 
housekeeping gene was employed as internal control. qRT-PCR was performed in 
a total volume of 50 fi\ containing: 1 x TaqMan Universal PCR Master Mix, no 
AmpErase UNG, and 2/A of cDNA using TaqMan assay on an Abi 7900HT 
Sequencing Detection System (Applied Biosystems, Paisley, UK). Transcript 
primers and probe sets employed were: GAPDH: primers FW: CTTTGTCAAGCT 
CATTTCCTGG, RW: TCTTGCTCAGTGTCCTTGC, probe CACCCTGTTGCTG 
TAGCCGTATTCA; BECN1: primers FW: GTACCGACTTGTTCCCTATGG, RW: 
ACACAGTCCAGAAAAGCTACC, probe: CCCC AG AAC AGTATAACG G CAAC 
TCC; CST3: primers FW: CTG ACTGTCCTTTCC ATG ACC , RW: TCCTTCTA 
GACTCAGCCCTTAG, probe: CCTTCCAGATCTACAGCGTGCCC; GABRA5: 
primers FW: GGGAATGGACAATGGAATGC, RW: TCTCATTGGTCTCGTCTTGTAC, 
probe: CATTTGCGAAAAGCCAAAGTGACCTGGA (Integrated DNA Technologies, 
Coralville, IA, USA). Real-time amplifications included 10 min at 95°C followed by 48 
cycles of 15 s at 95°C and 1 min at 60°C. For each transcript, relative expression 
levels were normalized against the GAPDH gene. The AACt method was used to 
compare relative fold expression differences. One-way ANOVA was performed to 
assess statistically significant expression differences among the analyzed 
transcripts. 
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